Short-interaction quantum measurement through an incoherent mediator 
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We propose a method of indirect measurements where a probe is able to read, in short interaction 
times, the quantum state of a remote system through an incoherent third party, hereafter called 
mediator. Probe and system can interact briefly with the mediator in an incoherent state but not 
directly among themselves and, nevertheless, the transfer of quantum information can be achieved 
with robustness. We exemplify our measurement scheme with a paradigmatic example of this 
tripartite problem: qubit-oscillator-qubit, and discuss different physical scenarios pointing out the 
associated advantages and limitations. 

PACS numbers: 03.65.Ta, 42.50.Ex, 42.50.Dv 



c ■ 

=5 

cr 



> 

(N 
(N 



On 
O 



X 



I. INTRODUCTION 

The transfer and extraction of information from quan- 
tum systems is an essential building block of quantum 
measurement theory This topic plays a key role 

in modern applications, like quantum computation and 
quantum information 0, quantum teleportation Q, and 
quantum cryptography \A\ . These applications have been 
driven by novel concepts in quantum information theory, 
as well as by successful measurement techniques achieved 
in different physical systems. Among them, we can men- 
tion electron-shelving qubit readout in trapped ions [f|, 
the quantum nondemolition measurement of field pho- 
ton numbers and the measurement of Rydberg atom lev- 
els in microwave cavity QED [6j. Recently, significant 
advances in the technology of quantum circuits Q have 
allowed efficient qubit readout techniques (8l-[l0| . From a 
theoretical point of view, there has been also recent in- 
terest in novel quantum measurement techniques, some 
of them aiming at improving experimental reach. We 
could mention weak measurements [ll| in photonics [l2j 
and circuit QED (13J, and short-interaction time mea- 
surements in the context of trapped ions 

EMI, 

cavity 

QED [HQ]!, and circuit Q ED GHH- 

When a direct interaction of a quantum system with a 
probe is available, the system information can be trans- 
ferred to the probe degrees of freedom. In this case, 
exchange of quantum information requires system-probe 
entanglement, but it does not warrant full reconstruc- 
tion of the system quantum state or even the proper en- 
coding of a searched system observable. In other cases 
we have to rely on indirect measurements, where a third 
party system mediates the communication between the 
system and the probe. In this sense, we are interested in 
the following question: is it possible that the probe mea- 
sures fully the quantum state of a remote system in short 
interaction times if only their independent coupling to 
an intermediate incoherent mediator is allowed? In this 



work, we answer positively to this question and discuss 
the physical frame behind it. In Sec. II, we introudce 
our model, while Sec. Ill describes a method of measur- 
ing indirectly the quantum state of a remote system by 
encoding information in the short-time behavior. Both 
unitary and dissipative cases are studied. In Sec. IV, we 
propose a possible realistic implementation in trapped 
ion setups. Finally, we discuss in Sec. V possible errors 
via numerical analysis. 



II. THE MODEL 

We consider the following model: a probe denoted by 
P will carry out the measurement of system S through 
a mediator A. We may also consider the mediator in 
contact with a reservoir, as depicted in Fig. [TJ The asso- 
ciated total Hamiltonian reads 



H = H P + Hs + H A + H P _ A + Hs-; 



(1) 



where Hp, H$, and Ha, represent the free energies of the 
probe, system, and mediator, respectively, while the cor- 
responding interactions arc denoted by Hp- a and H$-a- 
No direct probe-system interaction is allowed, that is, the 
effective coupling betwen probe and system can only hap- 
pen via higher-order terms of the unitary evolution. 

Let us consider a qubit(probe)-oscillator(mediator)- 
qubit (system) setup such that the initial tripartite state 
is decoupled, p(0) = p p (0)(g)pA(0)(S>Ps(0). The unknown 
system state is described by the density operator 

p s (0) = Pll |1) <1| + p 12 | 1) (2| + p 21 |2) <1| + p 22 12) <2| , (2) 

which we want to measure with a probe in state p p (0) = 
Pss I S> (g I + Pgo I g) (e | + Peg | e) (g | + Peo | e) (e | . The free energy 
terms are Hp = -^-a^, Ha = hu} a a^a, and H$ — ^^a z s , 
while the interactions are of the Jaynes-Cummings type, 



Hp- 



P-A 



hgj 



[a+a 



H S -a = hg s (a+a + a s a 1 ). 
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III. SHORT-INTERACTION MEASUREMENTS g p g s [(a+a- 



A. Unitary case 

As is the case in different experimental scenarios, we 
will consider the measurement of the probe observable 
9 p = |c)(c|, P c (r) = Tr [p(r)|e)(e|], as a function of the 
dimensionless interaction time r = g p t. From Eq. (TT|), we 
derive 



and write 



dP c (r) 
dr 



dP e (r) 
dr 



ihg p 



■\e)(e\,H P - A ]) 



(3) 



- (a+a - a v a*) 



(4) 



For r = 0, and choosing the initial probe state p(0) 
where = (|g) + e*+\e))/y/2, we get 



dP e (T) 



dr 



(5) 



T=0 



This shows that any mediator quadrature, X^ = (a^e 1 ^ + 
ae~ l< ^)/2, can be measured via the first derivative of 
P c (r) at vanishing interaction time [l4], [H, 113, 13 • ^ ue 
to the absence of direct probe-system coupling, it is clear 
from Eq. (j4|) that the probe is unable to record informa- 
tion about the system in a first-order expansion of P c (t). 

From Eq. (!]), we can derive further the second time 
derivative of P e (r) at r = 0, 



d 2 P e (r) 



dr 2 



1 



T = 



(ihg p ) 



[\e)(e\,H P _ A ],H 



(.6) 



After using [|e)(e|, Hp-a] = % P (ffpO — o>a^), the r.h.s. 
of Eq. ((6]) contains five remaining commutators associ- 
ated with each term of H . The four nonvanishing ele- 
ments, discarding [g p (a p a — <J p a*), ^fcr z ] = 0, read 



[g P (cr+a - a p a + ), ^cr^] = -g p u p {a+a + a p a)), (7) 
[g p {cr+a - cr~a t ),a; a a t a] = g p ui a (a+a + a p a^), (8) 



[g P (cr+a - a p a)),g p (o+a + a p a f )] = 2g p l (cr+<T p + a z n), 

(9) 




FIG. 1. (Color online) Sketch of the model. P, A, S, and R, 
are the probe, mediator, system, and reservoir, respectively. 



a p a 



),{<j+a + a s a))]=g p g s (aj 



(10) 

These terms are of second order in an expansion of P c (t). 
We observe that the first two terms interfere destructively 
and cancel each other in resonance, uj p = uj a . The third 
term will be mediator independent when (a p ) = 0, that 
is when p gg = p oc . The fourth term is the only one that 
could map, in principle, system information onto probe 
degrees of freedom. Specifically, it associates system co- 
herences with probe coherences, given that p gc = p* g =/= 0. 
It is clear, then, that a good choice of an initial probe 
state is p p (0) = \+ $) (+ <j>\- Note that in most quantum 
optical setups, where the probe is a two-level atom, the 
fidelity in the generation of this state is very high. Never- 
theless, we will estimate possible errors in the last section 
of this work. Finally, it can be shown that 



d 2 P c (r) 
dr 2 



■1 



T=0 



9p 
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(Pi 2 e l0 + P2ie-^) 



(11) 

where S = u p — u a is the detuning between the probe and 
the mediator. If the field is incoherent, p — '^2,P n \n){n\, 
e.g. a thermal state with (X^) = 0, or under the resonant 
condition S — 0, we are left with 



d 2 P c (i 
dr 2 



= -1 



2 9 P 



(pne** + pne-**). (12) 



This result shows that, though the probe and the system 
are only indirectly connected via an mediator, incoherent 
or not, it is possible to encode the coherences of the sys- 
tem in the short-time behavior of the second derivative 
of the probe's P c (t). A proper choice of phase 4> in the 
probe state |+$) will allow us to read the real and imag- 
inary parts of the system coherences. We remark that 
an incoherent mediator does not prevent the system to 
implement a fast transfer of information to the probe, as 
can be seen in Eqs. (TTTT) and (fT2|) . 

From previous considerations, we cannot gain informa- 
tion about the system populations due to the structure 
of the interactions Hp- a and Ha-s- Note that these 
interaction terms only couple, in first-order, off-diagonal 
elements of each qubit and the mediator degrees of free- 
dom. Given that the expression in Eq. (|12[) is related to 
a qubit-mediator-qubit second-order process, the system 
coherences can be, in principle, mapped onto the second 
derivative of the probe population. Nevertheless, we can 
measure the system populations p\\ and P21 provided we 
move from the quasiresonant regime to the off-resonant 
limit, where the qubits are dispersively coupled to the 
mediator. For the sake of simplicity, we will consider 
here that both qubits have the same transition frequency, 
uj p = lo s , and are coupled to the mediator with the same 
coupling strength g = g p = g s . In this case, we can write 
a second-order effective Hamiltonian from Eq. ([1} in the 
dispersive regime, 6 ^> g^J (n), as 



H, 



cir : 



hg 2 



[{a z p + a* s )<Ja + (a+ + a+){a~ + a;)] (13) 
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Remark that in this second-order Hamiltonian, as ex- 
pected, direct coupling terms between the qubits are ef- 
fectively generated. In this case, following similar steps 
as the ones used to obtain Eq. (TT2"]) , we can derive 



d 2 P c (r) 
dr 2 



T = 



3 ( 

§2 tP22 



Pll)- 



(14) 



It is worth mentioning that the result of Eq. ([14]) . similar 
to Eq. pip , can also work with any particular mediator 
state. In fact, this insensitivity could be used to test 
whether the couplings are quasiresonant or dispersive. 
From Eqs. (fT"2|) and (fl4|) . we can then recover the full 
density matrix describing the system, see Eq. @, which 
is encoded in the short-time behavior of the probe pop- 
ulation. In the last section of this article, we will study 
these results with analytical and numerical tools. 



B. Dissipative case 

We will study now tougher scenarios where the pre- 
vious results may fail to work. Let us consider then a 
single-mode quantum harmonic oscillator as the media- 
tor with a thermal bath acting on it [2l| . The associated 
master equation in the Markov approximation reads 



1 71b r t t 

P = -r[H, p\ + "f—[2a l pa - aa*p - 
in 2 

+ — (nt, + l)[2apa T — a'ap - 



- paa' 
pa) a] 



(15) 



where n& is the mean number of bath thermal excitations, 
and we will denote as n a the mean number of thermal ex- 
citations associated with the initial mediator state. We 
will consider the effects of the reservoir on our proposed 
readout of the system coherences. Surprisingly, we find 
that the presence of the reservoir does not affect the sec- 
ond derivative of P e {j) 1 preserving the results displayed 
in Eqs. ([12"]) and (fl4l . It can be proven that the ear- 
liest bath contribution appears only in the third-order 
derivatives, 



d 3 P c (r) 
dr 3 



1 



{my 



Tr p(r) 



\e)(e\,H P _ A ],H],H 



1 



Tr\£p(T)[[\e)(e\,H P . A ],H]\, (16) 



where the terms in the r.h.s. account for the unitary and 
nonunitary corrections to the second-order derivatives. 
If the initial probe state is a|g) + /3|e), the correction to 
Eq. (fTTj) , for instance, reads 



2 )(n a - rib). 



(17) 



It is noteworthy to mention that this contribution of the 
thermal reservoir becomes null whenever |a| = |/?| or 
n a = fib- The first condition can be easily satisfied when 
preparing the initial probe state, as in though we 



will discuss it in more detail below. The second condition 
is related to the lack of initial thermalization between the 
mediator and the bath. The fact that the reservoir affects 
Eqs. (TT^I) and (fbfj) only in higher orders is a signature of 
robustness of our proposed measurement scheme. 

To make a final test of the introduced concepts, we will 
consider now the case of both qubits interacting with a 
bath, that is the mediator consists of a continuum of 
modes. This model can be described with the Hamilto- 
nian 



2 



2 <y z P + hJ2^kblb k 

k 

+Kj29k[(v7 + Vp)b\ + (cr+ + a+)b k ] (18) 



We consider, for the sake of simplicity, similar qubit 
transition frequencies, w = uj s = uj p , and similar cou- 
pling between each mode and the qubits. This model 
contains some features of the unitary model described by 
Eq. (TT|), only that now the single mode has been replaced 
by many. It should be clear that both qubits are not 
coupled directly but they will get effectively entangled 
via the mediator modes in higher orders of the evolution. 
However, it is also expected that these modes play the 
role of a dissipative bath, making difficult the exchange 
of quantum information between the qubits. As a further 
scope, we would want to remark that going beyond the 
lowest-order perturbative regime of Eq. (fT8l) , a wealth of 
known effects appear in quantum information [22| . After 
tracing the bath degrees of freedom and using the Markov 
approximation, we are left with the master equation [23j j , 

-Op + cr ^)( cr p + < j s)p-p(< 7 p + cr ^)( cr p + O] » 

(19) 

where H = %VJ + ^cr*. In this case, and with |+^) as 
the initial probe state, we can calculate 



d 2 P c (r) 
dr 2 



= j [(Pii-p 22 )+ 2(1 + p 12 e i * + P2ie-^)} , 
4 

where the dimensionless time is now t = Pt. We observe 
that the populations and coherences of the unknown sys- 
tem are revealed all together at short-interaction times. 
This result can be intuitively explained in the following 
manner. The reservoir model for both qubits assumes 
the presence of a continuum of modes with different fre- 
quencies, some of them closer to resonance and some dis- 
persively coupled to the qubits. Equation (l20l) reflects, 
in an averaged manner, the influence of both cases, and 
it can be considered as related to the results found in 
Eqs. (fT2")) and (fTl]) . In consequence, it is possible to read 
the quantum information of a system with a probe only 
by the fact of sharing a common reservoir, as expressed 
by Eqs. ([IT)]) and (12TH) . Though this result seems to be 
strong, we must bear into account the influence of more 



4 



realistic bath parameters. When r is large, as expected 
for a real bath, a given discrete step of the dimensionless 
time t = rt would require a reduction in the step of the 
lab interaction time t. In consequence, a large thermal 
bath acting on both qubits would make more difficult for 
the probe to extract relevant system information. 



IV. EXPERIMENTAL ISSUES 

The present ideas could be implemented and tested 
in different physical setups. Due to their advanced pre- 
cision in quantum measurements, we believe that cir- 
cuit QED and trapped ion experiments are good can- 
didates, and we give details for the latter. We con- 
sider two trapped and laser-manipulated ions coupled 
through the vibrational center-of-mass motion like, for 
example, 40 Ca + ions held in a linear ion trap [ID, [25| 
with center-of-mass oscillation frequency v. In this case, 
each qubit can be realized by using the electronic level 
|45i/2(to = —1/2)) = \S) as the ground state and meta- 
stable level \3D 5 / 2 (m = —1/2)) = \D) as the excited 
state [24| . as shown in Fig.[3J with a magnetic field lifting 
the degeneracy of the Zeeman levels. In order to imple- 
ment our measurement protocol, we consider monochro- 
matic laser pulses individually addressing a single ion of 
the ion crystal [26] with a frequency cj close to the qubit 
transition frequency ojq. After applying an optical RWA, 
the Hamiltonian describing the ion-laser interaction reads 

HL = ^ {\D)(S\e"'^^ t+ae ^e^~^ +U.c}j . 

(21) 

Here, f2 is the coupling strength, 6 = uj — ujq the detun- 
ing, and r\ = k(h/2mi>) 1 / 2 is the Lamb-Dicke parameter 
defined as the square root of the ratio between the recoil 
energy and the ion center of mass quantum vibration en- 
ergy. Operators a, describe annihilation and creation 
of vibrational excitations of the center-of-mass mode. 



a Jaynes-Cummings Hamiltonian 
irjMl 
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FIG. 2. Realization of a qubit in a single trapped 40 Ca + ion. 
The qubit states \S), \D) are encoded in a pair of Zeeman 
states interacting with a laser pulse of frequency ui and cou- 
pling strength Q. 



In the case where the laser is tuned to the red side- 
band resonance, c5 = — v, the interaction (l2"Tj) turns into 



H, 



{a)\S){D\e- 1 * -a\D)(S\e lv ), (22) 



when the vibrational RWA is carried out and the ion is 
assumed to be in the Lamb-Dicke regime (j]{n) -C 1). 
Here, the phase tp corresponds to the relative phase be- 
tween the optical field and the atomic polarization. For 
ip = 0, the Hamiltonians Hp-A or Hs-a of Eq. ([I]) are 
realized. Measurement of the probe ion-qubit as required 
by (fT2|) or (fT4|) is efficiently carried out via electron shelv- 
ing [5[. In an ion trap experiment, the coupling strength 
g = rjtl needs to be kept much smaller than the oscil- 
lation frequency v in order to avoid off-resonant excita- 
tion of the ion on the carrier transition. For this reason, 
P c (t) changes on a time scale that is slow compared with 
the time scale oc v on which off-resonant excitations oc- 
cur. Measurement of first and second derivatives of P e {r) 
have been already accomplished in the lab with high pre- 
cision [H El]. They were implemented by measuring 
P c {t) for different probe times and fitting the measured 
signal by a low-order polynomial. Measurement of P e {T~) 
for a fixed probe time r entails preparing and detecting 
the ions a large number of times to keep the quantum 
projection noise at an acceptable level. Detecting, for 
example, 1% changes of P c (t) requires more than 10 4 
experimental repetitions which could be realised within 
a few minutes in current experiments. 



V. DISCUSSION 

In this work, we deal with (infinitesimal) derivatives 
of an expansion of P e (r) under some initial conditions 
and approximations. From this point of view, we should 
be careful when comparing our theoretical results with 
measured discrete derivatives in the lab. These consid- 
erations are not very different from the usual ones in 
similar models, but special attention should be given to 
derivatives. For example, Eq. (jf 2[) is valid only if we are 
not allowed or do not want to observe the probe dynam- 
ics with a temporal resolution beyond the rotating-wave- 
approximation (RWA). In this case, the second derivative 
of Eq. (fT2"|) predicts a correct measured value in the lab 
if the minimal discrete time interval of the experimen- 
tal sampling is At > g p /uj p . Otherwise, we would have 
to recalculate the model beyond the RWA to make the 
correct predictions. Similar arguments follow for the sec- 
ond derivative of Eq. fT4")) . Here, it is not only the RWA 
discreteness that should be requested, it is also the adi- 
abaticity of the second-order effective Hamiltonian that 
imposes an even slower temporal step in the dispersive 
limit. In consequence, the results of Eqs. (TTTI) and (TT4")) 
cannot be compared when 8 is large. In this sense, we 
have carried out numerical simulations to study and test 
the dispersive case for a particular example. Both, the ab 
initio Hamiltonian of Eq. (p} and the dispersive Hamilto- 
nian of Eq. (| 13[) were used, under similar conditions, to 
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FIG. 3. (Color online) Plots of P e (r e ff) and its derivatives as 
a function of the dimensionless effective time, T e B — g p t/S, 
for a dispersive case where the system state is \ip s ) = 0. 1 1 1) + 
e i7r/ Vl - 0.1 2 |2>, the probe state is \ip p ) = (|g) + \e))/V2, 
the detuning is S — 30g p , and the mediator is in a thermal 
state with n a = 1. The fast oscillating blue lines correspond 
to a calculation done with the ab initio Hamiltonian of Eq. (jTJ 
and the green line to the dispersive Hamiltonian of Eq. (1131) . 

test the predictions of our proposed protocol. For the ab 
initio model, the blue lines of Fig. [3] show the fast oscil- 
latory evolution that cannot be observed in experiments 
due to limitations in the temporal resolution. For the ef- 
fective model, that of Eq. (jT5)) . the green lines show the 
expected measured curves and derivatives. Note that the 
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FIG. 4. (Color online) Plots of the relative erros in the mea- 
surement of the system coherence and population, e pi2 and 
e P22 , respectively, and the infidelity, 1 — F, of the measured 
density matrix as a function of the population inversion AP. 
The state of the system is \tp s ) = 0.3|1) + 71 - 0.3 2 |2) and 
that of the probe is \ip p ) — (|g) + e})/\/2 . 



inset of the last plot in Fig. [3] shows the measured value 
of the second derivative, 0.975, at t c s = 0. 

We have also tested the robustness of our proposed 
scheme when deriving the first key result in Eq. (|12[) . 
We have considered imprecisions in the required vanish- 
ing of the population inversion of the probe, AP = (er^) 
as displayed in Fig. |4j In this figure, we have plotted the 
relative errors in the measurement of the system coher- 
ence and population, e Pl2 = (pi2 — /°i2 P )/pi2 and e P22 = 
(/J22 — P22 P )// 9 22, respectively, as a function of AP. Fur- 
thermore, we have plotted the infidelity of the measured 
state, 1 — F, also as a function of AP. Here, the fidelity F 
was defined with the help of the Frobenius norm and in- 
ner product as F = Tr[p cx V]/ V Tr [P cxp P oxp W Tl W]- 
From the plots in Fig. |4j as a key feature, one can observe 
that an imprecision of 0.01 in the population inversion of 
the probe will lead to a fidelity much better than 99%, 
assuring the protocol robustness under AP fluctuations. 
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VI. CONCLUSIONS 



cussed key features of the model and some fidelity issues. 
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